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Transporting and v T in a supernova (SN) core involves several processes that 
have been neglected in traditional simulations. Based on a Monte Carlo study 
we find that the flavor-dependent spectral differences are much smaller than is 
often stated in the literature. A full-scale SN simulation using a Boltzmann solver 
and including all relevant neutrino reactions confirms these results. The flavor- 
dependent flux differences are largest during the initial accretion phase. 



1. Introduction 

A supernova (SN) core is essentially a blackbody neutrino source, but in 
detail the fluxes and spectra depend on the flavor. Up to very small de- 
tails v^, v T , and v T can be treated on an equal footing and will be 
collectively refered to as v^. Numerical simulations usually find a hierarchy 
(E Ve ) < (Ep e ) < (E Vfi ) and approximately equal luminosities. The spectral 
differences offer an opportunity to observe flavor oscillations as the source 
fluxes will get partially interchanged. For example, it may be possible to 
distinguish a normal from an inverted neutrino mass hierarchy 1,2,3 ' 4 . 

A full-scale numerical simulation by the Livermore group finds for the 
integrated signal (E Vc ) = 13, (E 9c ) = 16 and (E u ) — 23 MeV and almost 
perfect equipartition of the luminosities 5 , results that are representative for 
traditional numerical simulations. Sometimes extreme spectral hierarchies 
of up to (Ep e ) : (E v ) w 1 : 2 have been stated, but searching the literature 
we find no support for such claims by credible simulations 6 . 

Traditional numerical simulations treat the and v T transport some- 
what schematically because their exact fluxes and spectra may not be im- 
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portant for the explosion mechanism. When a number of missing reactions 
are included one finds that (Ep e ) and (E V/i ) are much more similar than 
had been thought. The remaining spectral and flux differences are probably 
large enough to observe oscillation effects in a high-statistics galactic SN 
signal, but the details are more subtle than had been assumed in the past. 

2. Mu- and Tau-Neutrino Transport 

The transport of v e and v e is dominated by v e n <-> pe~ and D e p <-> ne + , 
reactions that freeze out at the energy-dependent "neutrino sphere." The 
flux and spectrum is essentially determined by the temperature and geo- 
metric size of this emission region. Moreover, the neutron density is larger 
than that of protons so that the v e sphere is deeper than the v e sphere, 
explaining (E Ve ) < (E„ e ). 

For i/p, in contrast, the flux and spectra formation is a three-step 
process. The main opacity source is neutral-current nucleon scattering 
v^N — ► Nvp. Deep in the star thermal equilibrium is maintained by nu- 
cleon bremsstrahlung NN <-> NNv^v^, pair annihilation e~e + <-> v^v^ 
and v e v e <-> v^D^, and scattering on electrons v^er — > e^v^. The freeze- 
out sphere of the pair reactions defines the "number sphere," that of the 
energy-changing reactions the "energy sphere," and finally that of nucleon 
scattering the "transport sphere" beyond which neutrinos stream freely. 
Between the energy and transport spheres the neutrinos scatter without 
being absorbed or emitted and without much energy exchange, i.e. in this 
"scattering atmosphere" they propagate by diffusion. 

One may think that the spectrum is fixed by the medium temperature 
at the energy sphere so that (E^) < (E v ) because the energy-sphere is 
deeper and hotter than the v e sphere. However, the scattering atmosphere 
is more opaque to higher-energy neutrinos because the cross section scales as 
E%, biasing the escaping flux to lower energies. For typical conditions (E v } 
of the escaping flux is 50-60% of the value characteristic for the temperature 
at the energy sphere 7 . Therefore, the final {E„ e ) : (E Vfi ) ratio is the result 
of two large counter-acting effects, the large temperature difference between 
the energy sphere and the v e sphere on the one hand, and the energy- 
dependent "filter effect" of the scattering atmosphere on the other. 

Until recently all simulations simplified the treatment of trans- 
port in that energy-exchange was not permitted in i/7V-scattering, e~e+ 
annihilation was the only pair process, and z^e-scattering was the only 
energy-exchange process. However, it has been recognized for some time 
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that nucleon recoils are important for energy exchange 6,7 ' 8 , that nucleon 
bremsstrahlung is an important pair process 6,7 ' 8,9 ' 10 ' 11 ' 12 , and more re- 
cently that v e v e — > v^v^ is far more important than e~e + — > v^v^ as a 
source reaction 6 ' 13 . 

We have performed a detailed assessment of the relevance of the new 
reactions on the basis of a Monte Carlo study 6 . To illustrate the results 
we use a hydrodynamically self-consistent accretion-phase model and show 
in Fig. 1 (left panel) the flux spectrum when using the traditional in- 
put physics (bottom curve). Then we add nucleon bremsstrahlung that 
increases the flux without changing much the average energy. Next we 
switch on nucleon recoils that depletes the spectrum's high-energy tail with- 
out changing much the overall particle flux. Finally, switching on v e v e 
annihilation increases the flux without affecting the spectrum much. The 
compound effect of the new processes is not overly dramatic, but so large 
that all of them should be included in serious full-scale simulations. 

In the right panel of Fig. 1 we compare for the same model the flux 
spectra of u e , D e and v^, the latter including all reactions. In this example 
(Ep e } almost exactly equals (E„ ) , but the fluxes differ by almost a factor 
of 2. This is reverse to the usual assumption of a pronounced hierarchy of 
average energies and nearly exact equipartition of the luminosities. 

We have studied a variety of stellar background models, some of them 
self-consistent hydrostatic models, others power-law profiles of density and 
temperature. For realistic cases we never find extreme spectral hierarchies, 
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Figure 1. Neutrino fluxes for an accretion-phase model. Left panel, curves from bot- 
tom to top: Flux of with traditional neutrino interaction channels, then adding nu- 
cleon bremsstrahlung, next adding nucleon recoils, and finally adding u e P e annihilation. 
Right panel: Fluxes for all flavors; the curve includes all reaction channels. 
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the differences between (E De ) and (E v ) typically being 0-20%. On the 
other hand, the fluxes can be rather different, especially during the accre- 
tion phase when the atmosphere is quite extended. The different neutrino 
spheres have then rather different geometric extensions, explaining large 
flux differences. Later during the Kelvin-Helmholtz cooling phase the star 
is very compact so that any geometry effect of the radiating surfaces is 
small. Moreover, the relevant regions are then neutron rich so that the 
transport physics of v e and will become similar. Therefore, during the 
late phases one expects very similar v e and fluxes, while, of course, the 
v e flux and spectrum remain unaffected by our arguments. 

3. Spectral Characteristics 

To characterize the neutrino fluxes one naturally uses some global parame- 
ters such as the particle flux, the luminosity (energy flux) , and the average 
energy (E). In order to characterize the spectral shape in greater detail 
one may also invoke higher energy moments {E n }. One measure frequently 
given from numerical simulations is E lms — (E 3 )/(E) because of its rel- 
evance for calculating average neutrino-nucleon interaction rates. 

Sometimes a global analytic fit to the spectra is also useful. Frequently 
one approximates the flux spectra by a nominal Fermi-Dirac function 



with a temperature T and a degeneracy parameter 77. This approximation 
allows one to fit the overall luminosity and two energy moments, typically 
chosen to be (E) and (E 2 ). However, the Fermi-Dirac fit is not more 
natural than other possibilities; certainly the low- and high-energy tails of 
the spectra are not especially well represented by this fit. 

We find that the Monte Carlo spectra are approximated over a broader 
range of energies by a simpler functional form that we call "alpha fit," 



E 2 



(1) 



1 + exp(-rj + E/T) 



f(E) oc E a exp [-(a + 1) E/E] . 



(2) 



For any value of a we have (E) = E, a Maxwell-Boltzmann spectrum 
corresponds to a = 2. The numerical spectra show values of a — 2.5-5, i.e. 
they are "pinched." 
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4. A Full-Scale Simulation 

In the Garching SN code 14 we have now implemented all relevant neu- 
trino interaction rates, including nucleon bremsstrahlung, neutrino pair 
processes, weak magnetism, and nucleon recoils. Our treatment of neutrino- 
nucleon interactions includes nuclear correlation effects. The transport part 
of this code is based on a Boltzmann solver. The neutrino-radiation hy- 
drodynamics program enables us to perform spherically symmetric as well 
as multi-dimensional simulations, thus allowing us to take into account the 
effects of convection. 

To explore the time-dependent properties and long-time evolution of 
the neutrino signal, we currently continue a state-of-the-art hydrodynamic 
calculation of a SN into the Kelvin-Helmholtz neutrino cooling phase of the 
forming neutron star. The progenitor model is a 15 M Q star with a 1.28 Mq 
iron core (Model sl5s7b2 from S. Woosley; personal communication). The 
period from shock formation to 480 ms after bounce was evolved in two 
dimensions. The subsequent evolution of the model is simulated in spher- 
ical symmetry. At 150 ms the explosion sets in, driven by neutrino energy 
deposition and aided by very strong convective activity in the neutrino- 
heating region behind the shock (Fig. 2). Note that a small modification of 
the Boltzmann transport was necessary to allow the explosion to happen 15 . 
Unmanipulated full-scale models with an accurate treatment of the micro- 
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Figure 2. Trajectories of the mass shells in the core of an exploding 15 Mq star. The 
explosion occurs about 150ms after shock formation, developing a bifurcation ("bubble") 
between the mass that follows the outgoing shock and the mass that settles on the nascent 
neutron star. Also indicated are the positions of the neutrino spheres of u e , P e and u T . 
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physics currently do not obtain explosions 16 . Details of this run will be 
documented elsewhere; at the time of this writing the CPU expensive cal- 
culation is still on the computer. Here we show in Fig. 3 a preview of the 
main characteristics of the neutrino signal up to 750 ms post bounce. 

The neutrino signal agrees with what is expected for the standard 
delayed-explosion scenario. In particular, it clearly shows the prompt v e 
burst and a broad shoulder in all fluxes during the accretion phase that 
ends at 200 ms when the explosion has taken off. The average neutrino 
energies follow the usual hierarchy and they increase with time due to the 
contraction of the star. We also show the alpha parameter from a global 
fit according to Eq. (2). During the accretion phase the flux is least 
pinched, at late times the a values of all flavors converge near 2.5. 

These results agree with and nicely illustrate our previous Monte Carlo 
findings in that the spectral hierarchy between v e and f M is rather mild and 
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Figure 3. Neutrino fluxes and spectral properties for the full-scale simulation described 
in the text. The hydrodynamic bounce and shock formation occur at t = (cf. Fig. 2). 
The right upper plot gives the spectral fit parameter E, the right lower one a. Note that 
the discontinuity in the latter at t Ri 480 ms is caused by mapping the model from two 
dimensions to one. 
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in that the average energies converge at late times. Conversely, the particle 
fluxes differ by almost a factor of 2 during the accretion phase, but cross 
over shortly after the explosion. At 750 ms the differences between the 
fluxes continue to increase, an asymptotic value has not yet been reached. 

5. Conclusions 

Traditional numerical SN simulations had two weaknesses regarding the 
flavor-dependent neutrino fluxes and spectra. First, the interaction between 
Vft.r and Dfi.r and the stellar medium was schematic, neglecting a number 
of important processes. Second, a Boltzmann solver for neutrino transport 
could not be coupled self-consistently with the hydrodynamic evolution. 

We have performed a systematic Monte Carlo study on various stellar 
background models and the first SN simulation that includes all relevant 
interaction rates and a Boltzmann solver. While the usual relationship 
between the v e and v e fluxes and spectra remains essentially unchanged, 
the Up spectrum is much more similar to that of v e , especially during the 
Kclvin-Hclmholtz cooling phase. Differences of the average energies are in 
the range 0-20%, with 10% being a typical number. During the accretion 
phase the particle flux is smaller than that of v e by up to a factor of 2, 
but later the particle fluxes cross over. 

Our findings imply that observing neutrino oscillation effects in a SN 
signal is a more subtle problem than had been thought previously, but by 
no means impossible. However, when exploring the physics potential of 
a future galactic SN one should not rely on the notion of an exact flavor 
equipartition of the luminosities or the extreme spectral differences that 
have sometimes been stated in the literature. 
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